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Abstract 

Since laws of physics exist in nature, their possible relationship to terrestial growth is introduced. 
By considering the human body as a dynamic system of variable mass (and volume), growing under 
a gravity field, it is shown how natural laws may influence the vertical growth of humans. This 
approach makes sense because the non-linear percentile curves of different aspects of human physical 
growth from childhood to adolescence can be described in relation to physics laws independently of 
gender and nationality. Analytical relations for the dependence of stature, measured mass (weight), 
growth velocity (and their mix as the body mass index) on age are deduced with a set of common 
statistical parameters which could relate environmental, genetics and metabolism and different 
aspects of physical growth on earth. A relationship to the monotone smoothing using functional 
data analysis to estimate growth curves and its derivatives is established. A preliminary discussion 
is also presented on horizontal growth in an essentially weightless environment [i.e., aquatic) with 
a connection to the Laird-Gompertz formula for growth. 
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I. INTRODUCTION 



The raw (and theoretically smoothed) statistical data of different aspects of human phys- 
ical growth from birth to adolescence have been collected for decades as a function of gender 
and nationality. From a practical point of view, these measurements are useful to monitor 
health care and to highlight secular trends on the fat intake (obesity) by a particular pop- 
ulation . Child growth records are also of special interest to Governments to enforce 
national health policies 3|. 

In particular, there are smoothed growth charts for boys and girls such as those main- 
tained by the USA National Center of Health Statistics (NCHS) ^\ and the unsmoothed 
data for thousands of Japanese infants, children and adolescents (Hiroshima Growth Study 
Sample) Typically, the charts consist of a set of non-linear percentile curves displaying 
the dependence on age of height h{t), "weight" w{t) and combinations of them such as the 
weight-for- height and body mass index (BMI, i.e., ratio B{t) = w/h?). According to the 
CDC growth charts [4], the 85*'* percentile of BMI for children is considered the overweight 
threshold, and the 95*'* percentile is the obesity threshold. 

In constructing smoothed statistical growth curves a variety of asymptotic mathematical 
models have been tested to fit results for a population in retrospective Jg^. Different trial 
functions are used because the data for h and w do not increase monotonically with age and 
because the use of weight-height methods (such as the BMI) affect directly the development 
of curve smoothing. For example, modeling of human anthropometric data has been done 



in the contexts of the lambda-mu-sigma ( 



^MS) model {?!, the model functions of triple 
9|, [lO| and Jolicoeur et al. curves 11[ . Another 



fine test model used is the infancy-childhood-puberty (ICP) model 12|. This model breaks 
down growth mathematically into different (exponential, quadratic and logistic) functions 
using different curve fitting procedures out of the growth data. Predictions of changes in 
height have also been carried out according to an empirical Bayesian approach 1^. The 
selection of a definitive parametric approach for the inclusion or exclusion of some empirical 
data points and to set a criteria of estimation is a topic of active discussions because of its 
relevance to growth, development and aging in living organisms js]. 

For open systems like those in which there are influx of mass, one can apply the con- 
servation of linear momentum and energy methods of general physics to study their growth 
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dynamics. Examples in which momentum is gain ed from, or lost to, the surroundings include 
rockets fl^. [isl and the falling of a snow ball 16 1, respectively. Since the human body is also 



a system of fluctuating mass (and volume), under the influences of the acceleration of gravity 
g and food consumption, it is not unreasonable to consider it as a dynamic physics system. 
One may then in principle make an attempt to derive a relationship based on physics laws 
and observed "physical" growth. 

Since laws of physics exist in nature, their possible relationship to human physical growth 
is introduced in this work. It is shown how physics seems to influence the physical growth 
of humans. Using Newton's second law, analytical relations are deduced for the time de- 
pendence of h, w, the growth velocity v = dh/dt which depend on common parameters that 
could relate other aspects of this phenomenon such as environmental conditions and genet- 
ics and natural processes like metabolism, energy supplied by food, etc. Within a simple 
physics-based framework, estimates of all observed human growth statistical data including 
BMI on age are easily performed and are shown to fit the observed data which makes the 
present model reasonable. A preliminary discussion on horizontal growth in an essentially 
weightless environment is also presented with a connection to the Laird-Gompertz formula 
for aquatic growth. 

II. SIMPLE PHYSICS CHARACTERIZATION 

Let us consider the human body as a system of variable mass where mass is the 

amount of matter present in the body measured using balance. In Fig.l the approximated 
center of mass of an upright body as seen from a particular reference frame at different ages 
is illustrated. In scientific usage mass is an intrinsic property of matter and weight is a 
force that results from the action of gravity on matter. In everyday usage, however, weight 
and mass are used interchangeably. The reported "weight" in the smoothed growth charts 
for humans in 4, ISj corresponds to measured matter or mass in kilograms and not to the 
force of weight "w{t) = m{t)g" measured in Newtons. In the following such "weight" data 
is referred as measured mass for the sake of correctness. Changes of mass and momentum 
will be assumed to be continuous during the growing process through years from childhood 
to maturity. 

To derive the equation of motion of the body system whose mass is not constant, let 
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us use Newton's second law to define the external force F on the human body. As shown 
explicitly in Appendix A, this is expressed as (see, e.g., [l4| ) 

dw ^ dm 

m — = Fext + ^rel -J- , (1) 

dt dt 

where dm/dt > is the rate at which mass is gain [i.e., equivalent to particles entering 
an open system) and Vrei{t) is the relative velocity of the gained mass with respect to 
the body moving with a growth velocity of magnitude v{t) = |v(t)|. The above term 
(f)[t) = \rei {dm/dt) corresponds to the momentum flux being tranferred into the body by 
the added mass. This is interpreted as the force exerted on the body by the mass that joins 
it [l^ and F^xt is the external force of gravity acting on the variable mass system. The 
surface of the earth is taken as the zero level of gravitational potential energy. 

In applying the Newton law for the vertical growth of humans, it corresponds to the 
analysis of an extended system of particles rather than to the analysis of a single particle 
of mass mi and velocity Vj at time t. In this case one can only consider averaged values 
for all the physical quantities -such as the external force, due to the complexity of the body 
system. Hence, if n{hi) is the number of objects at various heights hi, of a total of n objects 
distributed among i heigths, then the mean average value of any given function / for those 
objects is in general 

1 ^ ^ n(h ) ^ 

where pi is interpreted as the probability that a fraction of the objects n{hi)/n has the 
height hi. So as the sum of all fractions of the whole system is unity (/ = 1), the probability 
distribution function is also normalized such that Y^\=i pi = 1. It can be seen that pi directly 
affect the outcome of averaged results since it weights the values of the given function / at 
each height i. The above mean value equation is used throughout the theory of statistical 
physics. 

Let us consider next the number of n objects distributed among i heights as being A^ 
sets or groups of particles of added mass mi each at time t. Then the mean value for the 
external force of Eq.([l]) in one longitudinal direction becomes 



^ „ ^ i dvi dmi\ . 
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Statistically stationary states in which all the probabilities p are time independent Vi will 
be only considered. 

To an observer at rest on the earth, the accumulating particles in the body at a given time 
appear to move at a velocity, say fa(t), proportional to the growth velocity v{t) of the center 
of mass. For simplicity let us consider an homogenous system so that such a proportionality 
is set equal to one, i.e. Va{t) ~ v{t), which corresponds to the maximum velocity that the 
added mass alone can achieve during the grow. This is so because the human growth in 
stature is longitudinal until reaching a maximum height h{oo) with respect to earth's surface 
(and a maximum mass m(oo) under healthy conditions) at adolescence (see Fig.l). On the 
other hand, to an observer on the center of mass of the body, the added mass appears to 
have a vertical motion with a velocity of magnitude v{t) + Va{t) ~ 2v{t) in the opposite 
direction of the growth with respect to the center of mass displacements. Therefore, for the 
present purposes, the relative velocity Vrei(^) of the gained mass appears to move at twice the 
velocity of growth as seen from within the body system, namely \vrei\ = 2|v|. Hence Eq.(j3]) 
for the mean value of the net external force on the body as a whole is the superposition of 
contributions 

Fe^t) =< F,,Ut) >= Ep. m{t) ^-2J:p, v,{t) . (4) 

i=i i=i 

with p the normalized probability distribution function of N objects or set of particles of 
added mass rrii each at time t. 

The force of gravity Fg due to the earth mass acting thoughout the human body must be 
also taken into account. Similarly to suspended chains systems this force is not only 
proportional to the total weight of the body at an given height h{t), but also to the total 
force due to the amount of mass at rest lying on the surface of the earth. This extra term 
can also be understood by the fact that the volume (besides the mass) of the human body 
also varies. It accounts for an additional momentum that the system acquires as a result of 
changing volume (in one-dimension vertical) with years -see Appendix A. In this case, the 
force of gravity acting on the N objects of added mass rrii is then given by 
^ r h it) 1 

where Xi{t) = mi{t) / hi{oo) represents mass per unit height. The weight at rest term w{oo) = 
m(oo)g =< mj(oo) > g = Y^f=i PiTni{oo)g is deduced according to the boundary condition 
Fg ^ Q as t ^ oo. 



III. GENERAL BIOMETRIC FORMULAE 



In order to derive analytical relations for the dependence of stature, measured mass 
(weight), growth velocity (and their mix as the body mass index) on age, let us proceed by 
considering Newton's law of motion in relation to the two distinct mean forces acting on the 
system, i.e. the net external force of Eq.(jl]) and the gravity force of Eq.(j5]). This means to 
set Fg = Fgxt, or more explicitly 

f , , dh^it) (dhi{t)\ (dmAt)\ hAt) , , \ 

with Vi{t) = dhi{t)/dt. After some little algebra, and using the well-known trigonometric 
relation: sech'^ + tanh^ = 1, it is straightforward to check that this second order differential 
equation has the general solutions 

^ ft-t\ 
hit) =<h,{t)>=Y,Pihi{.^) tank { -\ , (7) 

i=l \ Ti J 

such that 

h.{^) = \gTl , (8) 

with tj, and time lags for i sets of particles. The mean value of Xj is deduced by considering 
the limit t — > oo, which leads to h{oo) =< hi(oo) >= Y^iLi Pihii^o) = g < rf > /2. In this 
way, the mean time lag < Tj > is interpreted as the time at which the growth velocity becomes 
zero and the maximun body height is achieved with respect to the earth's surface h{{)) = 
(recall that x{t) = x(0) + v{t) * t + |a(t) * from Newton's mechanics for constant mass 
systems). On the other hand, tj are the time lags at which the single height contributions 
hi{t) of the objects or set of particles become zero. 
Furthermore, 

v{t) = ^=<^>=j:p.^sech^('-^) , (9) 

which for t — > oo, it gives f (oo) ^ as expected (after completion of the physical growth), 
and 

^ ft-ti\ 
m{t) =< mi{t) >= ^ piiTLiioo) tank i ) . (10) 



Ti 



i=l 

It follows inmediately that physics theory predicts that the height and mass functions reach 
maximum values (plateau) on age due to the presence of the hyperbolic tangent functions. 



The mass per unit height Ai(oo) = mj(oo)//ij(oo) ^ corresponds to the slope of mj(t) 
versus hiif). For systems of variable mass, the center of mass velocity vit) is not the same 
as that for the N sets of particles in the system f In turn, the moving particles may not 
have the same velocities for each state i. 



IV. PHYSICS ESTIMATES OF HUMAN GROWTH DATA 



To this point, the amount of growth parameters in the present physics based biometric 
formulae seems to be large. However, these are not only comparable with those in the alter- 
native asymptotic mathematical models tested for the modeling of human anthropometric 
data (see, e.g., 0, Q, [sl, 1^, Q, H, Q, Q), but their number can still be reduced -at least 
for small sets of particles. 

It is easy to verify by a plot that for y > 0, the monotonicity of the hyperbolic tangent 
implies the relation (iTJ 



tanh y < 1 < 



tanh {xy) 



(11) 



tanh X 

Therefore, using this inequality it is possible to group together model parameters to set lower 
bounds to all functions in this work. In fact, by setting y = r/rj ^ l/cj > (independent 
of t) and X {t — t^jr, such that = e^r, one obtains 



tanh [xy] = tanh 



t-ti 



T 



> tanh ( — ) tanh 



t-t, 



T 



Hence Eqs.(I7]) and ffTOl) can be approximated as 



N 



h{t) >J2^i ^^^^ 



and 



where, Vi, it is defined 



1=1 



N 



t-ti 



m 



(t) > ^ Bi tanh 



i=l 



T 



t-ti 



tanh y tanh x . (12) 



(13) 



Ai = pihi{oo) tanh ( —- \ = pihi(oo) tanh ^ 



Bi = Xi{oo)Ai 



(14) 



(15) 



.Cj/ \Ti 

with A is mass per unit height as before. Since the Ti parameters are related to hi via 
the Newtown equation in Eq.(l8]), then this leads to the positive dimesionless quantities 

ei 



{1 / t) y 2hi{oo) / g . Furthermore, the momentum flux tranferred to the body b y t he 
added becomes (pi{t) = Vi {dnii/dt) = {Bi/Ai) Vi{ty, similarly to nonrigid systems [15|. 



Taking the derivative of Eq. (fT3!) it follows that 

1 ^ 

v{t)>-J2Aisece(- — -) . (16) 



which is consistent with Eq.dH]) in the lower bound since {1/Ti)sech'^{xy) > 
(1/r) tank y sech^ x. 

In this way, the unkown growth parameters hi{oo), mi{oo) and pi can reduce to Ai and 
Bi, which will be referred to as "biological parameters". Together with the time lag r and 
the peak po^it.ons of each tank function, given by i^the^e pa.a.,ete.= allow to describe 
non-linear percentile curves (or alternatively z-scores [6]) as shown next. 

Selected smoothed empirical percentiles of measured mass (weight)- and height-for-age 



from the boys data based on medical literature ^\ are shown in Figs 2 and 3 by dots, 
respectively. The theoretical curves (full lines) are obtained using Eqs. (fT4l) and (fT3l) assuming 
= 6 sets of particles with the biological parameters Ai/ Bi of Eq.l fTSi) illustrated in Fig.4. 
All parameter values can easily be derived using the fit command of the free GNUPlot data 



plotting software. Real (smoothed and unsmoothed) data curves 



monotonically with age independent of gender and race 



0, 0, y, y: 



or boys a nd g irls increase 
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12 



13|. In the 



inserts, unsmoothed 50% percentile empirical points for w{t) and h{t) of Chinese girls ages 



6-18 years 18| are shown for comparison to display observed characteristic trends. 

In order to show that there is not need for guessing some mathematical functions a priori, 
let us fit the smoothed curves rather than the raw data. A similar behaviour to the reported 
smoothed curves is obtained from the Newton dynamics of systems with variable m. Within 
measurements error, the variation of estimates of growth charts is close to all smoothed 
empirical data in the figures (usually measured to the nearest 0.1cm [19]). Both the height- 
for-age and measured mass-for-age tend to a maximum. According to these results, it can 
be argued that for each plotted point, the variable mass and height (and growth velocity as 
seen below) of humans can be correlated in relation to physics laws over many years. 

At the origin, the boundary conditions /i(0) = and m(0) = in Eqs.( fT3l) and f |T^ lead 
to the relations between the biological parameters 

Ai tank (^^^ - ^ ^^^^ (7) ' ^1 ^"^^ (7) - ^ ^^^^ ( r ) ' ^^^^ 

where ti < and t2, - ■ ■ ,tN > which leads to reduce further the number of parameter 
values. In fact, as t ^ 00, tank and B-^^ ^0 for finite N sets of particles (of 



added matter m.^^). 

In humans growth is a natural process characterized by periods of growth spurts during 
early infancy and again at puberty. In Figs. 5 and 6 illustrated by dots are the growth velocity 
curves and the measured mass (weight)-for-height curves at different percentiles, respectively 
^. The full lines in these figures correspond to results obtained from Eq. fll6l) and Eqs. fll4p - 
(fT3|) . using parameters as before. As shown in Fig. 5, the growth velocity oscillates down 
as a function of time to vanish in the adolescence. The system acceleration dv/dt changes 
sign over the years due to the tanh term in the derivative of Eq. (fT6!) . On the other hand, 
both the real and the mass- for- height model curves of Fig. 6 increase monotonically. The S*'* 
percentile in this case is considered as an underweight threshold [4] . Such weight-for-height 
standards are used in the study of nutrition status and mortality risk in a community |3|. 

Further evidence for a relation between physics laws and human growth is given by 
looking at the BMI on age plotted in Fig. 7 for selected percentiles. On equal grounds, the 
plotted dots are evaluations using the measured mass and heigth on age using the smoothed 
charts of Figs. 2 and 3 [S] and the full line curves from m{t)/h{tY =< rni{t) > / < hi{t) >^ 
with the biological parameters in Fig. 4. It can be seen that the present estimates of BMI 
are reasonably close to the reported smoothed data. BMI fluctuations increase with the 
error of the data for h"^. This is reflected in BMI calculations for younger ages at the higher 
percentiles. Theoretical BMI curves at the lower percentiles display a minima with increasing 
age from 3 to about 6 years to then rise steadily with age as observed in reality. Within 
the physics approach, this minimium is a consequence of the growth spurts at early infancy 
depicted in Fig. 5. Increasing in the sums of Eqs. fll4l) and (JT3l) allows a better minima fit 
at higher percentiles of the BMI-for-age curves. 



V. REMARKS 



The present study is not based on arbitrary mathematical trial functions to fit the hu- 
man physical growth data. All general trends for the observed non-linear percentile curves 
displaying the distribution of different aspects of human physical growth from childhood to 
adolescence have been understood and described in relation to physics laws that are daily 
perceived by living on the surface of the earth. 

As shown here theoretically (see Appendix A for details on the derivation of the general 
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Newton approach in question) , the underlying physics of Newtonian mechanics is the bases 
of the derivation of all equations. The coefficients at each stage of the derivations are the 
consequence of dealing with the differential equation ([T]) applied to a very dense data set 
based on human growth. This also include the fundamental function, tanh{x) appearing 
throughout the resulting relations for the body height in Eq.([7]) , mass in Eq. fllOp and BMI 
thereafter. Therefore tanh{x) here is not made arbitrarily to fit a wide range of data from 
a set of adjustable, weighted coefficients. On the other hand, the peaked growth velocity of 
Eq.Q depends on the different hyperbolic function sech{x). 

Alternative models to estimate growth curves and their derivates, such as the growth 
velocity, are based on the use of a non-parametric monotone smoothing process of the form 



The regression coefficients /5o,i set the origin and range required to fit /i data including the 
velocity curves v{t) on age to account for the pubertal growth spurt. The fi(t) — > v{t) 
function is constructed using the functional data object W{u) which does not approximate 
the measured data directly, but rather after some exponential transformations (to make 
fi positive) and integration up to t (to make /i increasing or monotonic). A relationship 
between the present physics based approach and the monotone smoothing can be established 
as shown explicitly in Appendix B. 

The tendency to predict and describe growth by normal curves often obscures individual 
growth spurt phases that are consistently present but asynchronous across the population 
as a whole ^|. To this end, one can deduce that since the present growth theory is based 
on the law of gravity which influences individuals, then either a single contribution, or the 
superposition of contributions, on the growth curves are driven by the same general patterns 
described as derived by Newton's law in Eq.([T]). 

When dealing with an extended system such as the human body, Newton's second law 
states that the net external force applied to the center of mass equals the center of mass 
acceleration times the mass, namely Fg^t = Ma(cm). According to Eq.(jl]) one can then 
write Fext =< rriiVi (dvi/dhi) — 2vf {drrii/dhi) >=< mj(oo)aj(cm) >. Using the average 
deflnition in Eq.(l2]), e.g.,m =< rrii >, in conjuction with the deflniton of mass-per-unit height 
Aj(oo) = mj(oo)//ij(oo) and Eqs. (fT3|) - (fT5|) it follows that hi(t) = {hi{oo) / mi{oo))mi(t) and 
Vi{t) = hi{oo){dmi{t) / dt) . Hence hi{oo)ai(cm) = hi {dvi/dt) —2 v^. The latter corresponds to 




(18) 
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a variable velocity f , namely L x a™^^) = v'^ + x {dv/dt), and that of a uniform rod of height 
L' = L/3 lifted with constant velocity v' = v, namely L' x = f'^. This superposition 



the difference between the center of mass acceleration of a uniform rod of height L lifted with 

, rod 
''{cm) 

''(cm) 

may reflect the fact that the center of mass of the human body -whose mass distribution 
can also change by pelvic rotation and displacement, knee flexion, foot mechanisms- has not 
a fixed location on time 2l|]. There are alternative models based on Newton mechanics to 



study mass redistribution throughout the body as a function of age and applied forces (see, 
l^l)- -'^y Ei'Pplyiiig the equation of motion for the human body with variable mass as 
discussed here, it could be helpful to understand the extent to which empirical models by 
inverse dynamics to derive the torque (moment) correspond to prediction as well as their 
affect on human movement. 

The primarily interest in human physical growth and its relationship to physics laws is 
to enhance the fit quality to observed data and to understand the underlying phenomena. 
There is overwhelming biological evidence that continuous growth processes on earth proceed 
by multiplicative increments, therefore the body size of plants, animals and humans follows 
a lognormal rather a normal distribution (with additive increments) 23|. Thus one should 

n 

not expect human height distributions to be normal but lognormal [24]. The claim here is 
not modified by either class of pi distributions -relating the biological parameters in Eq. (fT5!l . 
In principle, one may associate these parameters to changes in kinetic energy of the center 
of mass energy, via the integral of Fext{t) in 

Eq.dH) with respect to = x'^/2L - x'^/2L'. 

The relationship between natural laws and human physical growth as introduced in this 
work may be also suitable to understand the biometrics of other living species growing 



vertically such as Primates 



251 ] (and perhaps leaving trees and plants too). However, at 



this point, it is also worth to ask if the same physics laws discussed here can also be applied 
in the case of horizontally growing terrestial organisms, i.e., those organisms growing in an 
essentially weightless environment (aquatic organisms as white-side Dolphins 26|). 

So the question is: Can horizontal growth be also described within the present physics 
approach in the "absence" of gravity (g 0) as contrary to vertical growth where gravity 
influences directly? 

Let us attempt to reply this crucial question and see how the present model may provide a 
unique description for both cases. Consider the general Eq. (jlj) for Fext again and set the force 
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of gravity equal to Fg = 0. Alternatively, from Eq.dS]) this implies hi(t) hi{oo) = const. 
which means to compensate out the weight w{t) with the volume dependent forces. Therefore 

, ^ dvAt) , ^ dmAt) , . 

J2p,m,it)^~2j2p,v,{t)—^ = . (19) 

i=i i=i 



One solution satisfies 



or 



Vi{t) rriiit) 
which, by integration, implies the relation 



ram ^ = 2.,«) . (20) 



dvi{t) _ 2 dmi{t) ^21) 



-.W-^=7m?(t) , (22) 

or 

v{t) =< Vi{t) >= ^ < ml{t) > , (23) 

with 7 a constant for a given aquatic specie and for the long x-axis of a horizontal structure. 

On the other hand, let us consider (and expand in Taylor series) the well-known Laird- 
Gompertz empirical formula for fish growth [26|, |27J . That is 

^ 1 + a (1 - e-°*) 
1 



1 + a 1 



gat 



where x{t) is the lenght at age t (at zero gravity) and (a, a) are specific rates of the expo- 
nential growth. It follows that the growth velocity v{t) = dx{t)/dt satisfies 

vit) = j^^ = —[il + a)xo-xit)r . (25) 
(1 -|- ar)^ axQ 

Using the Chebyshev's sum inequality J2^=i^kbk > (1/''^) Sfc=i «fc Z]fc=i ^fc, such that 
ttk = bk = PiTTii = prrii, then comparison of Eqs.(|22l) and (1251) results in < mi{t) >^ < 
N pT.i=i pmiitf = Np< mfit) > and 

1/2 

axo^ 



rait)<(^] |(l + a)x„-xW| . (26) 
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which for values x{t) > (1 + a)xo, it imphes an approximated hnear relation between mass 
and length for horizontal growth (the mass-for-height for vertical growth under the influence 
of gravity is plotted in Fig. 6 for comparison) . ^ 

To this end, collected data reported in [26] for white-side Dolphins in the central north 
pacific ocean suggest that the postnatal dentine thinkness of Dolphins (used as an index 
of age in these animals) increases linearly with the total body length. It is reasonable to 
then justify the linear proportion between the aquatic species mass and length of the type 
in Eq. fl26|) . as a consequence of the application of physical laws (c./., Eq.(jl])), in conjuction 
with the Laird- Gompertz growth function (c./., Eq. (1241) ) . Relations ( 12^ and ( l26i) also allows 
to study the BMI for aquatic organisms and Eq. (1231) to predict their growth velocity curves 
in terms of the square of their mass (c./., Eq. fl2^ ). 

Therefore, the absence of gravitational forces during growth in secondarily aquatic 
tetrapods (SAT) -backboned, limbed, air-breathing animals with ancestors that evolved in 
a terrestrial environment where gravity was a factor in growth and motion {e.g., Dolphins)- 
leads to mass and body length trajectories that, in principle, could be replicated by the 
present general physics model. Further multidisciplinary research along these lines is needed 
to shed new light on the fundamentals of growth in all classes of terrestrial organisms using 
one model based on Newton's second law as discussed here. 
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APPENDIX A 



Let us follow Refs. Il5| to demonstrate Newton's law in Eq.([T]) for open systems 
where the mass varies with time and is not constant. 

The net external force acting on the system can be approximated as 

_ rfP AP _ P/ - Pi 

~ ^ ~ At ' ^^^^ 

where Pf^i denotes the final and initial system momentum, respectively. As an illustrative 
example, let us consider a mass m moving with velocity v that ejects a mass Am during a 
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finite time interval At. One can then write 



Pi = mv ; P/ = (m - Am)(v + Av) + Amu , (28) 

where u is the velocity of the center of mass of the ejected mass (not longer in the original 
system). The system mass is reduced to m — Am and the velocity of the center of the system 
is changed to v + Av. Therefore, replacing Eq. (!28l) into (1271) leads to 

^ Av ^ , . Am , . 

Fe.t =m — + [u-(v + Av)] ^ . (29) 

The quantity u — (v + Av) is just the relative velocity \rei of the ejected mass with respect 
to the main body. 

The change in mass with time, is intrinsically negative in this example. As At 0, then 
the positive quantity Am/ At can be replaced by —dm/dt. Hence, 

dv dm . . 

m- = F.., + v..— . (30) 

The term Vrei{dm/dt) is the rate at which momentum is being transferred into (or out of) 
the system by the mass that the system has ejected (or gained). It is interpreted as the 
force exerted on the system by the mass that leaves (or joins) the system. 

When the volume is allowed to vary, a term accounting for the additional momentum 
that the system acquires as a result of its changing volume (distict from the momentum flux 



(f){t) = \rei (dm/dt)) needs to be added to Eq. flSU]) |28| . In the model for BMI this is given 
by the volume term in Eq.Q. 



APPENDIX B 



The computation of W{u) is difficult and not unique. Its basis expansion needs to be 
estimated iteratively. Some tips on monotone smoothing as given in {2^ suggest to use 
n-order B-spline basis functions for W{u). Within this functional data analysis (FDA) it is 
necessary to know a-priori the right behaviour of ^(t) in order to calibrate the smoothing of 
growth data with the smallest possible number of knots. 

A relationship between the present physics based approach and the monotone smoothing 
using FDA is given next. Let us expand Eq. f|T6|) for the growth velocity in Taylor series for 



14 



sech{x), such that 

1 ^ / 



24 



1 ^ 



1 ^ 



^ -J2Aexp{-x{ty} ; (Ixd < vr/2 < 00) . (31) 
i=i 

In the above 

xW-(^) , (32) 

and 4t/i-order terms in the expansion of sech{x) and exp{x) ^ 1 — (or ahernatively terms 
of the Gth-oidei since, to a good approximation, (2/3)x^ ~ x^/2! are neglected. 

In the continuous hmit {i.e., J2i ~^ I du), the discretized growth velocity function in 
Eq. (!3T!) -which varies in terms of the i index from 1 to A^, it becomes 



v(t) ^ — / exp{~u'^} du = — < / exp{— M^} du+ exp{—u'^} du > 

T J{t-ti)/T T [Jt/T Jo J 

(33) 

The integral on the right of the sum corresponds to the well-known error function erf[{t — 
tis[)/T\ for finite t <t]^ and large r, with t]^ the expected value and r the standard deviation. 
For simplicity an uniform growth due to the added mass objects and approximated ti — 
and y4j — > A is assumed. 

Therefore a comparison of Eqs. (JTSl) and (!33ll imply these associations 



T 

2 



W{u) -u' . (34) 

Therefore, a link to FDA and the present approach such that v{t) > follows by setting 
y4 < and — = W{u) = Y17=o ^i,pi'^)Pi where the sum gives the B-spline function of 
degree p with Pi the contact points and u the knot vector. The function W{u) in Eq. (!33|) 
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leads to the continuous normal (or Gauss) distribution with zero mean and variance of one. 
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FIG. 1: Representation of human physical growth as a dynamical physics system of variable mass 
(and volume). 



18 




I • 1 • ' • 1 • 1 • ' • 1 ' 1 • ' • 1 

24 48 72 96 120 144 168 192 216 240 

Months 



FIG. 2: Measured mass (weight)-for-age for boys at selected percentiles. Real data (dots) using 



smoothed charts in 



J] and theoretical results (full lines) from Eq. ()14p . Unsmoothed 50% percentile 



empirical points for girls (of a different nationality) 



181 ] are shown in the insert for comparison. 
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FIG. 3: Height-for-age at percentiles as in Fig. 2. Real data (dots) using smoothed charts in 4| 
and theoretical results (full lines) from Eq. (ll3p . 
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FIG. 4: Inverse of the mass per unit height parameter Aj = 1, • • • , A^) of Eq. (ll5p used throughout 
the physics based calculations. Full lines are used to guide the eyes. 
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FIG. 5: Growth velocity curves on age at percentiles as in Fig. 2. Real data (dots) using smoothed 
charts in [4] and theoretical results (full lines) from Eq. ()16p . 
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FIG. 6: Measured mass (weight)-for-height at percentiles as in Fig.2. Real data (dots) using 
smoothed charts in P] and theoretical results (full lines) from Eqs. (jl4p - (jl3|) . 
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G. 7: Body Mass Index at percentiles as in Fig.2. Real data (dots) using smoothed charts in 
J] and theoretical results (full lines). 
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